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Abstract. An integrable self-adjoint 7-point scheme on the triangular lattice and an integrable self-adjoint scheme 
on the honeycomb lattice are studied using the sublattice approach. The star-triangle relation between these systems 
is introduced, and the Darboux transformations for both linear problems from the Moutard transformation of the 
B-(Moutard) quadrilateral lattice are obtained. A geometric interpretation of the Laplace transformations of the 
self-adjoint 7-point scheme is given and the corresponding novel integrable discrete 3D system is constructed. 
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> ■ 1. Introduction 

O 

1.1. Integrable linear systems on triangular and honeycomb grids. In this paper we study integrable linear 
problems on planar graphs built of regular polygons (of the same type) different from squares: the self-adjoint linear 
problem on the star of the regular triangular lattice, and the self-adjoint linear problem on the honeycomb lattice. 
We call a linear problem integrable, if it possesses a certain number of relevant mathematical properties, including: 
the existence of i) discrete symmetries (Laplace and Darboux type transformations), ii) continuous symmetries 
(nonlinear evolutions), iii) dressing procedures (9-dressing, finite gap constructions) enabling one to construct 
large classes of analytic solutions. 

To understand the integrability of systems on the triangular and honeycomb lattices is a challenging subject. 
Exactly solvable models of statistical mechanics on such grids are rather well studied [B], but there are not so 
many papers discussing this problem from the point of view of integrable difference equations (see, however, 

da EH n 13 H [a El 1311 HE] ) • 

^ I The self-adjoint linear problem on the triangular lattice was introduced into the theory of integrability in |32, 34] 
\ ' as a discrete analog of the two dimensional elliptic Schrodinger operator. This discretization was distinguished, 
C — ■ . among other operators having the same continuous limit, by the existence of a decomposition into a sum of a 
' multiplication operator and a factorizable one. Such decomposition, or "extended factorization", leads to trans- 
, formations, called in the literature Laplace transformations [32 [M|- The transformations preserve the form of the 
' operator but does not involve any parameters. A more general transformation of the Darboux type (with functional 
parameters) for the self-adjoint linear problem on the triangular lattice was introduced and studied in [28j . 

Integrable circle patterns on the plane with the honeycomb combinatorics were studied in [71 [SI H] not only 
from the point of view of integrable systems, but also as potentially important objects of the discrete holomorphic 
function theory. In [22] special Laplace transformation operators on the triangular lattice were studied, within that 
context, as analogs of the d and d operators. In [7] circle packings exhibiting the so called multi-ratio property 
^ ' were investigated, while the patterns with constant intersection angles were studied in [8]. A special reduction of 
the packing with constant angles leading to the discrete Painleve and Riccati equations was the subject of [S]. 

1.2. Main ideas and results. The theory of integrable difference (or discrete) equations on the grids is 
much more developed. It turns out that many integrable systems can be obtained by reductions of the 
multidimensional quadrilateral lattice [TH [T51 [T71 [1^] (multidimensional lattices with all elementary quadrilaterals 
planar), called also discrete conjugate net [37]. One of the goals of this paper is to show that the integrable systems 
on triangular and honeycomb grids can be incorporated into the multidimensional quadrilateral lattice theory. 
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Figure 1. The diagonal staircase section as the bipartite quasiregular rhombic tihng (on the left), 
and the corresponding triangular and honeycomb lattices (on the right) 

By imposing an additional linear constraint on the quadrilateral lattice, one obtains the so called B-quadrilateral 
lattice [14j, which is an important object in the paper. On the algebraic level, such a lattice is described by the 
linear system of discrete Moutard equations [TTI l3T|, the compatibility condition of which gives the nonlinear 
Miwa equations [27], called also the discrete BKP equations (discrete Kadomtsev-Petviashvili equations of type 
B). Recently, integrability properties of the two-dimensional version of this lattice were used to explain integrability 
of the self-adjoint 5-point scheme on the star of the 71? lattice [U]. The approach presented there was based on 
the idea of deriving from integrable equations on a lattice (and from their integrability features) novel integrable 
equations on a sublattice of the original lattice, together with their integrability features. In particular, i) the 
Darboux transformation for the self-adjoint 5-point scheme given in }28j was rederived in |21| from the Moutard 
transformation _31j of the discrete Moutard equation, ii) the finite-gap theory for the self-adjoint 5-point scheme 
was obtained as direct consequence of the finite gap theory (constructed also in [5T]) for the discrete Moutard 
equations in Z^. 

This is also the main idea behind the present paper. We derive the integrability properties of the self-adjoint 
linear problems on the triangular and honeycomb lattices from the corresponding properties of the three dimensional 
B-quadrilateral lattice. In particular, we obtain in this way the Darboux-type transformations for both linear 
problems, and we show that the Laplace transformations of [32j are nothing but a relation between two subsequent 
triangular sublattices of the B-quadrilateral lattice. We remark that the approach we use here is very close, in 
spirit, to the one used in where properties of the hexagonal circle patterns with constant intersection angles 

were derived from properties of an integrable system on the T? lattice, and to that used in [lOj . where discrete 
holomorphic functions on quad-graphs were constructed from discrete holomorphic functions on 1^ . Our results 
can be also interpreted from the point of view of the discrete Cauchy-Riemann (discrete Moutard) equation on 
quad-graphs I26j. with the graph being the quasiregular rhombic tiling, as visualized on Figure [T] However, we 
would like to stress the multidimensional origin of our constructions and their geometric meaning. 

The layout of the paper is as follows. In Section O after construction of the triangular and honeycomb grids 
as sublattices of the diagonal staircase section of the 1? lattice (see Figure [ij, we derive the self-adjoint 7-point 
scheme and the honeycomb linear problem from the discrete Moutard system. The sublattice approach provides 
the clear geometric meaning to the star-triangle relation between both linear problems. In Section [3] we present the 
corresponding interpretation of the Laplace transformation of the self-adjoint 7-point scheme [32j as the transition 
between subsequent triangular sublattices of the B-quadrilateral lattice. We introduce a 3D fully discrete novel 
integrable system describing such a transition. Finally, in Section [4] we obtain the Darboux transformations of the 
self-adjoint 7-point scheme [28j and the honeycomb linear problem from the Moutard transformation. In the rest 
of this introductory section we recall necessary informations on the discrete BKP equation, its linear problem (the 
discrete Moutard system), and the corresponding Moutard transformation. 

1.3. The discrete Moutard system (the B-quadrilateral lattice). Consider the system of discrete Moutard 
equations (the discrete BKP hnear system) [TTl [31] 

(1.1) V^, -i> = r X (^, - ^,), 1 < z < J < Af, 
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where ip : ^ N, and V is a linear space over a field K, and : ^ K, 1 < i < j < iV, are some 
functions; here and in all the paper by a subscript we denote the shift in the corresponding discrete variable, i.e. 
/i(ni, . . . ,ni, . . . ,njv) = /(^i, . . . ,ni + 1, . . .,71^)- 

Remark 1.1. In JUIJ it was shown that the discrete Moutard system p.ip characterizes algebraically (up to a gauge 
transformation) special quadrilateral lattices x = [V'] : P(V) in the projectivization of V, which satisfy the 

following additional local linear constraint: any point x of the lattice and its neighbours Xij, Xik and Xjk, for i, j, k 
distinct, are coplanar. 

The coefficients are not arbitrary functions. Compatibility of the Moutard system leads to the following set 
of nonlinear equations [31j 

(1.2) i + /f (r-r'=) = /fr^■ = /^■r^ i,j,k distinct, 

where, formally, we put = . 

If 6 : ^ K is a scalar solution of the system then the solution ijj : — > V of the Moutard transformation 
equations [U] 

(1.3) 4,,^iJ = U^-^^). l<^<N, 
satisfies a new discrete Moutard system with the fields 

(1.4) f''=f''j^- 

We will need also the fact that 9^^ satisfies the same Moutard system as the tranformed wave function ip. 

Remark 1.2. The discrete Moutard transformation provides the algebraic part of the B-reduction of the fundamental 
transformation of the quadrilateral lattice 14J. This transformation acts within the class of B-quadrilatcral lattices 
and, apart from the usual property of the fundamental transformation stating that the corresponding points x, Xi, 
X, Xi of the two quadrilateral lattices are coplanar, we have the following additional linear constraint: the points 
X, Xij, Xi and ij, for i ^ j, are coplanar as well. 

2. Linear problems on the triangular and honeycomb lattices 

The goal of this section is to present the derivation of the discrete Laplace equations on the triangular and 
honeycomb lattices from the system of Moutard equations on Z^ . This is done in full detail to prepare the setting 
for the new results in Sections[3]and|4l The derivation of the Laplace equations from the discrete Moutard (Cauchy- 
Riemann) equations on planar bipartite quad-graphs is well known [261 19]. The Z^ origin of some integrable systems 
on planar graphs was used in [51 [5] , where properties of the hexagonal circle patterns with constant intersection 
angles were derived from properties of an integrable system on the Z^ lattice, and in [lOj , where discrete holomorphic 
functions on quad-graphs were constructed from discrete holomorphic functions on Z^. 

2.1. The staircase section. Define the bipartition of the Z"^ lattice depending if the sum of the coordinates is 
even (black points) or odd (white points). Consider the subset 

V{T) = {m(l, 1, 0) + n(0, 1,1): m, n e Z} = {(ni, na, "2) G : ?^l - ?i2 + "3 = 0}, 

of the black point lattice, which can be considered as the intersection of the Z"^ lattice with the diagonal plane. 
Intersecting, instead, the Z'^ lattice with two "nearests" parallel planes we obtain the subsets 

V{T+) = {(1,0,0) + to(1,1,0) -|-n(0,l,l) : m,neZ} ^ {(^i, 712,712) G Z^ : 711 - 712 -|- rig - 1 = 0}, 

y(T_) = {(0, 1, 0) + 7?i(l, 1, 0) + 7i(0, 1, 1) : 771, 71 e Z} = {(ni, 712, "2) ^ : 711 - 712 + 773 + 1 = 0}, 

of the white point lattice. We call the union of the above sets, together with the corresponding edges and facets of 
the Z'^ lattice, the diagonal staircase section of Z^ (see Figure [T]). As a planar graph the diagonal staircase section 
is often refered to as quasiregular rhombic tiling (the dual of the kagome lattice). 

Points of V{T) are vertices of the triangular grid T whose set of edges E{T) consists of segments connecting the 
nearest neighbours 

EiT) - {[vi,V2] : Vi,V2e F(T), \vi - ^ V2}. 
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Figure 2. The triangular lattice, its main directions and the fields a, b, s on its edges 



Similarly we define the triangular grids T± , which will be used in Section [3l However we will need first the 
honeycomb grid H whose set of vertices is V{M.) = V{T+)U V{T-), and edges of which are the segments connecting 
points of V{T±) with their nearest l/(Tzp) neighbours, see Figure [T] 



We keep the bipartition of the honeycomb lattice into zt-points. The honeycomb tiling is the dual of the triangular 
tiling (in the standard crystallographic sense). Roughly speaking, this means that the centers of the facets of one 
lattice are the vertices of the second one (see Figure [T]). 

2.2. The self-adjont 7 point scheme and the sublattice approach. The goal of this Section is to derive, 
using the sublattice approach, the self-adjoint 7-point scheme from the system of the Moutard equations on 
grid (the three dimensional B-quadrilateral lattice). However, we first recall a way (see, for example [TU]) in which 
one can construct the Laplace equation for an arbitrary graph Q with given weight functions v : E{Q) — > K from 
non oriented edges to the field K. When tp : V{G) ^ V is a function on vertices of the graph, then it satisfies the 
corresponding Laplace equation if at each vertex Vq 



where the summation is over vertices adjacent to Vq. 

In the case of the triangular grid there are three families of edges, and therefore function i/ gives rise to three 
functions, say a, b and s, given on triangular lattice; see Figure [U where also the corresponding convention of 
labeling directions / and // of the lattice is given. The Laplace equation on such weighted graph reads 

(2.2) a(V'/ - -0) + «_/(?/;_/ - -0) + K'^ii - 0) + b^u{^-ii - -iA) + siiipi-u - ip) + sniip-iji - -0) = 0, 

where ip : V{T) ^ V is a field defined in the vertices of the grid with values in a vector space V. It is easy to see 
that given any solution : ^ V of the the Moutard linear system p.ip . then its restriction to the vertices V{T) 
of the triangular lattice satisfies the self-adjoint affine 7-point scheme (|2.2p with the edge fields defined as 



E{m) = {[v+,v 



(2.1) 





and the shifts on the triangular lattice given by 

(2.4) (/) = (12), and (//) = (23). 



Equation (|2.2p is the affine form of the general self-adjoint scheme on the triangular grid 
(2.5) A*/ + A-/*-/ + S*// + B_7/*_// + Si^i-ii + Sii'^-i.ii = F^, 



considered in [32]. In this paper we confine ourselves to consideration of the affine scheme (|2.2p (without loss of 
generality because, as we will see below, every scheme ()2.5p can be transformed to an affine form ()2.2p by means 
of gauge transformation). 
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Figure 3. The honeycomb lattice 



Corollary 2.1. Given a scalar solution a of the generic self-adjoint 7-point scheme (j2.5p . then tp = /a satisfies 
the affine self-adjoint 7-point scheme ()2.2p with the coefficients 

(2.6) a = Aaicr, b = Baua, s — Sa^ia^u . 

Conversely, given a generic function p : V{T) ^ K on the triangular grid and given a solution ip : V{T) W of 
the self-adjoint 7-point scheme (j2.2p . then the function ^ = ip/ p satisfies the scheme (j2.5p with the coefficients 

{a + a^i + b + b^n + si-n + s^iji)p^. 



(2.7) 



A 



apip, 



B = bpup, S = sp^jp^ii, 



F 



Remark 2.2. Since the Moutard equation on the grid can be integrated by the finite gap method [2X1, the above 
sublattice derivation of the 7-point scheme from the discrete Moutard equations imply that also the 7-point scheme 
on the triangular lattice can be integrated using the theory developed in [5TI [13] . 

Remark 2.3. Notice that the above derivation of the self-adjoint affine 7-point scheme can be generalized to arbitrary 
dimension as follows. Consider a solution ip : ^ V of the the A^-dimensional Moutard linear system ()f on 
the 2N elementary quadrilaterals 

(1, 2), . . . , ij,j + 1), . . . , (TV - f , N), (N, -f ), (-1, -2), . . . , {-N, 1). 

Add the first N schemes and substract the remaining ones (dividing by the Moutard coefficients, when necessary). 
The result is the following self-adjoint (2iV + f )-point scheme (on a star with 2N legs) on the black points: 



2.3. The linear problem on the honeycomb lattice. The discrete Moutard equations (jl.ip imply a linear 
relation between a point of the T_ lattice with its three nearest T+ neighbours: 

(2.9) /12(V-i - + /''(^3 - ^2) - 7^(^-123 - V^2) = 0, 

J2 



and a linear relation between a point of the lattice with its three nearest T_ neighbours: 

1 



(2.10) 



r"(^2-V'i) + /ii3(V' 



12-3 



^1 



713- (^-3- 
^-3 



0. 



These two relations, which give the geometric characterization of the B-quadrilateral lattice [14] , can be interpreted 
as the following linear self-adjoint system of equations on the honeycomb white lattice (see Figure 

(2.11) 1(0+ - r) + \{<t>ti u - r) + —(01/ - r ) = o. 



(2.12) 



1 



SIJI 

1 



br-, 



= 0. 
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for the fields and (j) , defined by 

(2.13) 0+=i/'i, r=V'2 



having used also (|2.3p and the triangular lattice shifts 

Again, like in the triangular lattice, one can consider a more general linear system on the honeycomb lattice 
than the above affine one (I2.1ip - (j2.12p . The result below is the analogue of Corollarv l2.1l for the honeycomb lattice. 

Corollary 2.4. Given a solution ((7^,a^) of the generic honeycomb system 

(2.14) ^$+ + ^$1,^^^+ 1^^^$+ =F+<i>-, 

(2.15) >" + ^'^^"- + i*=--^"^"' 

then ((/)+, (/)^), where cjy^ — $^/cr^, satisfies the affine honeycomb linear system (|2.1ip - (|2.12p with the coefficients 

A , B S 



(2.16) 



Conversely, given generic function {p^ , p ) : 1^(]H[) —> on the honeycomb grid, and given a solution {(f>'^,4) ) : 
^(H) —fWof the linear problem (|2.1ip - ()2.12|) . then the function ($+,$") = ((/)+/p+, ^^/p^) satisfies the system 
(P7Ti)l - (PTT^ with the coefficients 



P^P P-i,iiP P-iP-i.-ii 

\a b siji J \a bi^^u sj 

2.4. The duality between the linear problems on the triangular and the honeycomb lattices. As we 

have already remarked, the honeycomb and triangular tilings are dual to each other. In fact, the duality can be 
extended to the level of the linear problems. Both the honeycomb linear system (|2.1ip - (|2.12p and the self-adjoint 
affine 7-point system ()2.2p are consequences of the following " duality" equations between function on the T lattice 
and the pair (</>+ ,(/>") on the honeycomb lattice 

(2.17) = 

a 

(2.18) ^jj-^ = l{^- -cptjjj), 

(2.19) V//-V'/ = — (0|/-r)• 

s/,// 

Notice that the star-triangle equations (I2.17p - (l2.19p are just the Moutard equations on the diagonal staircase 
section. Equations ()2.17p - ()2.19p have a point of T_ as the star center, while equations (I2.17p . ()2.18P r -fr and 
(|2.19P - rr have a point of T+ as the star center (see Figure H]). This is the direct counterpart of the famous star- 
triangle relation widely used in the analysis of solvable models in statistical mechanics 0. As it was shown in 
pi] the (nonlinear) star-triangle transformation in electric networks, interpreted as the local Yang-Baxter equation, 
leads to Miwa's discrete BKP equation. The (linear) star-triangle relation (|2.17p - (|2.19p gives thus the linear problem 
for the nonlinear one (see also Remark l3.3p . 

Remark 2.5. The above star-triangle relation can be considerd as a particular case of the relation between the 
Laplace equations obtained from the discrete Cauchy-Riemann (discrete Moutard) equation on quad-graphs [9], 
where the graph is the quasiregular rhombic tiling. However, we would like to stress the three dimensional (not 
planar) origin of our construction and its geometric meaning. 
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Figure 4. The star-triangle duality relations between the triangular and the honeycomb tilings; 
the dashed lines represent the edges of the quasiregular rhombic tiling (the staircase diagonal 
section) 



3. The Laplace transformation of the triangular lattice 

The Laplace transformations of the self-adjoint 7-point scheme on the triangular lattice were defined in (32J and 
studied in |341 122j . The construction used in |32j consists in the factorization (up to a multiplication operator 
summand) of the second order operator (|2.5p into product of two first order operators. In this Section we will 
present the geometric meaning of the Laplace transformations. The basic idea is as follows: the hyperplanes 

= {(rti, rt2, «2) G I?^ ni — n2 + — £}, ^ £ Z, 

foliate the bipartite lattice into triangle lattices of a fixed (black for £ even, and white for £ odd) colour. The 
transition between two subsequent lattices of the same colour is the Laplace transformation. 

To develop algebraically this idea we will use formulas (|2.1ip - (|2.12p of the affine honeycomb linear system, which 
can be considered as the connection formulas between the affine self-adjoint 7-point schemes on T±i = T± lattices. 
To simplify our calculations we will work in the affine gauge induced by the Moutard linear problems (in [32 El] 
generic gauge was used). Notice that, in the (forward) transition from to T_|_, we have three natural possibilities 
(see Figure S]) 

This ambiguity and the corresponding one for the backward transition was observed on the algebraic level in [33]. 
In this article we choose the transition 0~ ^+ (the others can be obtained by simple shifts) as the forward 
Laplace transformation £+. 

The second equation ()2.12p of the honeycomb linear problem can be rewritten in the form 

(3.1) = £+0" = ^ (^bi-iisjcf)^ + abi^^ii(pZu + "-^i^l-ii 
where 

(3.2) r ^ a^js + b^jjs + a^jb^jj. 

Remark 3.1. The above transformation is the afhne version of the Laplace transformation P2 of |33]. 

Inserting such into the second equation of the linear system (j2.1ip - (|2.12p . we obtain that the field (/)~ satisfies 
the affine self-adjoint 7-point scheme with the coefficients 

/o ON - «ff 1.-^1 - ^ 

(3.3) a = , = , s = — . 

riJi riji r 
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Remark 3.2. The coefficients are given up to a common constant multiplier. The above choice is motivated by the 
sublattice approach. Indeed, equations ()2.13p and ()2.3p imply that, in the sublattice notation 

n- - — h~ ~ — - - f 1^ 

J2 J2 



We leave it for the Reader to check that equations p.3p are consequences of the discrete BKP equations 

Remark 3.3. Equations p.3p are directly related to the star-triangle transformations in electrical (or spring) net- 
works, see for example |24, between the conductivities of conductors (stiffness of the springs); compare also 
Section [231 



The above equation can be reversed, i.e., the coefficients a, 5, s can be expressed in terms of a , 6 , s as 

o-Zn s~ _______ 1 

(3.4) a = — , 5— — s = , where q =a b +as^jj + bsjjj — 



I-11 1-1 1-1 ^-n ' ' ""i^ii 

This allows to write down the forward Laplace transform of (j)^ in terms of coefficients (13. 3p of the scheme on T_ 

(3.5) 0+ = £+0" = ^— (bZuSjcl)^ + aZibZu<PZii + aZjSj (jyj . 

l-ii ^ ' ' 

Similarly, the backward transformation of 0+ 

(3.6) 0" = £"0+ = i (bsjji(j)+ + ab(l>jj + asiji4)+jjj'^ , 
where 

(3.7) q = asiji + bsjji + ab, 

inserted into the ffist equation of the linear problem gives the affine self-adjoint 7-point scheme with the coefficients 

(3.8) a+ = -^^, b+ ^ 



qi-ii q q-ii 

Remark 3.4. Again, the above coefficients can be obtained by the sublattice approach, i.e., 

Jl Jl 

and equations p.Sp are consequences of the nonlinear discrete BKP equations (|1.2p . 
The inversion formulas 

(3.9) a = T— -, 6=—-^, s = J^, where = a'tjb'tii + a'tjS^ + b'tjjS^ = , 

rjj rjj rZj q-ii 



allow to write down the backward Laplace transform of 0+ in terms of coefficients p.8p of the scheme on T+ 
(3.10) r =C-c^+ = \ (b+sijj4>+ + at J jjb+^jj + at J u4 ii^-i 11) ■ 

We have obtained therefore the sublattice derivation of the Laplace transform, which we formulate in terms of the 
proposition below. 

Proposition 3.5. Given a solution of the self-adjoint affine 7-point scheme with coefficients p.3p . then its 
Laplace transform (f)'^ = C'^<j)~ , given by formula p.Sp . satisfies a new self-adjoint affine 7-point scheme on T+ 
whose coefficients can be obtained from those ofT^ by inserting expressions p.4p for functions a, 6, s into equations 
(|3.8p . The inverse Laplace transform from 0+ to (j)^ — (jf^ is given by equation p.lOp . 



Corollary 3.6. Interpreting the Laplace transformation as a shift in a new discrete variable -fC £ Z, from equations 
(|3.4p and p.9p one can write down a nonlinear system relating coefficients of the subsequent 7-point equations 
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where 
(3.12) 



7^ - a^6^ 



^ -^1,11 







7^ 



Remark 3.7. As we have seen, the discrete (3D) BKP equations ()1.2p can be interpreted as the connection formulas 
(or the discrete time evolution) between subsequent self-adjoint problems on the triangular lattices T^, ^ G Z. This 
implies also that one can obtain the full Mutard system on the 1? lattice from the self-adjoint linear problem on 
the triangular lattice (or the honeycomb lattice). 



Remark 3.8. The novel, up to our knowledge, nonlinear system p.lip - (|3.12p is the sublattice version of the discrete 
BKP equations. It describes the transformation rule of the coefficients of self-adjoint 7-point schemes under the 
Laplace transformation. The corresponding transformation rule in the case of the discrete affine hyperbolic 4-point 
scheme is equivalent to Hirota's discrete Toda system [53], and was derived in [T^] (see also p3|)- geometric 
counterpart is given by the Laplace transforms of two dimensional discrete conjugate nets, introduced in pure 
geometric way by Sauer |36| . The Laplace transform of the hyperbolic 4-point operator (in general gauge) was also 
defined independently, using the factorization technique, in 32 . 

We would like to close this Section by presenting the decomposition of the affine self-adjoint 7-point scheme in the 
spirit of the papers [3H[22] (formulas below are the affine versions of the corresponding formulas in there). Denote 
by Ti and T// the standard translation operators along directions I and // of the triangular lattice. Combining 
both Laplace transform p.ip and p.6p we obtain equations 

^ (6/,„//S/ + abi^^iiTjj^ + asjTiTfj^) ^ {bsiji + abTu + asijiTjTjj^) 



- (bsiji + abTu + asuiTjTjj^) 



{bi-iisi + abi^^ijTjj^ + asjTiTjj^) 



i.e., "factorized" linear problems on and T_. Notice that the fields a, b, s are the coefficients of the self-adjoint 
7-point scheme ()2.2p of the intermediate black-point lattice Tq . 

Elaborating the above formulas, one can derive the following decomposition of the affine Laplace operator of the 
lattice (p^ (we leave the case of 0+ to the Reader) 



(3.13) 

in terms of the operator 
and its formal adjoint 
where 



[Q-{Q- 



Q- 
{Q- 



y Tii 



0, 





a b 




^1,11 


1-II 






aZn 




^1,11 



Remark 3.9. In [31] the relation between the coefficients p. 31) of the 7-point scheme on T_ and the corresponding 
coefficients of Q~ has been brought to the form 

VJ- = {x^f + {y-f + {z-f + a- + aZi + b^ + bZu + ^7 + 

Remark 3.10. In the special case x — y = z = 1 the operator Q and its adjoint Q* were viewed in |22j as the 
analogs of the operators d and d in the discrete holomorphic function theory. In the Dirac operator on more 
general bipartite planar graphs was defined and related to the dimer model of statistical mechanics. 



Remark 3.11. The Laplace transformation equations p.ip and (|3.6p are nothing but the linear problem (|2.1ip - (|2.12p 
on the honeycomb lattice. Conversely, the factorization (j3.13p implies the honeycomb linear problem. 
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Figure 5. Geometry of the sublattice approach to the Darboux transformation of the self-adjoint 
7-point scheme 



4. The Darboux transformations 

Our goal here is to derive the Darboux transformation of the self-adjoint 7-point scheme [55] and of the linear 
honeycomb problem using the sublattice approach. In doing that we will follow the corresponding ideas of |21j . 
where the Darboux transformation of the self-adjoint 5-point scheme |28l was rederived in such a way. 

4.1. The Darboux transformation of the linear problem on the triangular lattice. Before presenting 
the sublattice derivation of the Darboux transformation, let us first comment on the geometric content of the 
calculations below. The geometric definition of the B-quadrilateral lattice states that there exists a linear 
relation between ■0, -0-1-2, "0-1-3 and ■0_2-3; i.e., the corresponding points [0], [0_i_2], [0_i_3], [■0-2-3] of 
the projectivization of V are coplanar. Similarly, the geometric definition of the corresponding reduction of the 
fundamental transformation states that there exists a linear relation between ■0_i_3, 0_2-3, 0-3 and 0_i_2-3. 
Both corresponding planes intersect along the line passing through [■0_i_2] and [0_i„3], which means that all the 
six points belong to a three dimensional subspace. Therefore any five of the six points are linearly related (see 
Figure ini). This allows, in particular, for the elimination of 0_i_3, which does not belong to the allowed subset 
V{T) of the Z'^ lattice, yelding a linear relation between -0, ■0_i_2, '0_2-3, "0-3 and '0_i_2-3, the first equation of 
the Darboux transformation. 

Three discrete Moutard equations shifted in (— z — j) directions, lead to the linear equation 

(4.1) 7^— (V' - ^--1-2) - 7T5— (V' - V'-i-a) + Ta^CV' - ^-2-3) = 0. 

/-I-2 /-I-3 /-2-3 

li 6 is a scalar solution of the system then the Moutard transformations (|1.3p for i — 1,2, taken with 

appropriate shifts, and the discrete Moutard equation, shifted in the (—1 — 2 — 3)) direction, give 

^ — 2 " ]^ ~ ~ 

(4.2) (0-3 - V'-2-3) - -5 (V'-s - i/'-i-a) - ^7:; (V'-a - 0-1-2-3) = 0. 

-3 



-2-3 (^-1-3 ./ 1-2-3 



The elimination of "0-1-3 from equations ()4.1|) and (14. 2p leads to 



~ ]_ ~ ~ y"-— "1 — 3 ^ — 3 

(4.3) 1p-3 + ^ (0_3-0_l-2-3)--7 (0'_3-"0-2-3) = "0-712 (V' " ^"-1-2) - 723 ("0 " "0-2-3) 

^^-3 /_i-2-3 c'-2-3 /-1-2 /-2-3 

We multiply both sides by /l'i_2-3^-2-3/li-2//-'i-3, ^-nd we use (a part of) the nonlinear equations H 

(4.4) - /il-2/"l-2-3, 
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to transform equation (j4.3p into 

(4.6) 0_3^-3 - e-l-2-3V^-l-2-3 = -^zl^V^ + /ii_2-3^-2-3V'-l-2 + f /ii_2-3^-l-2 ) V'-2-3. 

J -2-3 \J-2-3 

Similar considerations, with equation 

(4.7) 7r^(V^-i - ^-1-3) - 7r^(V^-i - ^-1-2) - -7^ — (^-1 - ^^-1-2-3) = 
t^-1-3 /-1-2-3 

instead of equation (|4.2p . lead to 



(4.8) - 0_i_2-3V;-l-2-3 - TT^^V- - /li_2-3^^-l-2V'-2-3 " 7T2 /-i-2-3^-2-3 ) ^'-1- 

J-1-2 \j-l-2 



2- 



Introducing the field ip — 6'_i_2-3V'-i-2-3 and using the sublattice notation (|2.3p . the above system (|4.6p and 
is rewritten in the form 



(4.9) ^i-ip ^ -h^uB^ni' - sO-h^j-i + (6_//0 + sB^i)i'^u, 

(4.10) 'ipii-^^ a^iO^iip + se^ii^^ii- {a^iB + se^ii)ip^i. 

Because and 9^^ satisfy the same Moutard system, then they satisfy the same self-adjoint affine 7-point scheme. 
By Corollary 1 2.11 the function tp satisfies the self-adjoint affine 7-point scheme whose coefficients can be found from 
equations (|2.3p . with f^^ in place of supplemented by the (—1 — 2 — 3) shift and by the change (|2.6p . with 
respect to the above-mentioned gauge, with a = l/0_i_2-3- This implies that the coefficients read 

1 1 



6= - 



/1^2-3^-3^-l-2-3 n\-2-3(^-2-30-l-3 ' 
1 1 



/'!-2-3^-l^-l-2-3 /-i-2-3^-1-2^-1-3 ' 



£ — _ f 13 ^ — _ f 13 ^ 

^III - -1-1-3 n n - •'-l-3flfl ■ 
[7-1(7-3 [7W_1_3 

To rewrite them in the sublattice notation we notice that equation (14. 5p can be rewritten in the form 

1 r 



(4-11) .13 
where r is given by l|3.2p . Then equation (|4.4p and its analogue for /^^ give 

(4-12) TT2 = 723 = ~- 

J -1-2-3 ' J-1-2-3 ' 

Finally, equation (|4.ip . satisfied by 9, and equation (|4.1ip imply 

^-1-3 , 9a^ib^u 

fl3 

J -1-3 



which gives 



(4.13) = a_i9_i + b_u9-n 

9-uP 9-ip Bp 

where 

(4.15) p = s{a^i9^i + b^u9-ii) + a-ih^u9. 

Finally we can formulate the theorem, which was given (in the equivalent non-affinc form, see Remark l4. 21 below) 
in [55] directly on the 7-point level. 

Theorem 4.1. Given the general solution of the affine self-adjoint 7-point scheme (j2.2p . and given its particular 
solution 9, then equations (14. 9p and (|4.10p define the general solution tp of the affine self-adjoint 7-point scheme 
with coefficients given by formulas (j4.14p - (j4.15p . 
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Figure 6. Geometry of the sublattice approach to the Darboux transformation of the hnear 
problem on the honeycomb grid 



Remark 4.2. Equations (|4.9p and (|4.10p and the transformation formulas (|4.14p for the coefficients coincide with 
equations (32) and (34) in [28]. The identification ip — tp'O, where t/j' is the corresponding transform in [28], implies 
that tp' satisfies a self-adjoint 7-point scheme (an arbitrary gauge does not change the structure of the equation) 
which will not have, in general, the affine form. 

Remark 4.3. One could consider another (equivalent) choice of the Darboux transformation direction, for example 
V' —> ^I23'0i23- This would change slightly the Darboux transformation equations (|4.9p and (|4.10p and formulas 
([4T4l) - (l4l5l) . 

Remark 4.4. Notice that, for ip defined on T^(T), the sum of the ni, n2 and tlt, coordinates of ^-1-2-3 is odd. 
However, it is convenient not to considered '!/'_i_2-3 as a field defined in points of the white lattice. The proper 
interpretation would be to consider the Moutard transformation direction as a shift in the fourth dimension of the 
Z** lattice. Then the sum of all coordinates (including 724 = 1) of "0-1-2-3 is even. 

4.2. The Darboux transformation of the integrable honeycomb lattice. In this last part of the paper we 
first state the theorem on the Darboux transformation of the honeycomb linear problem, and then we show its 
sublattice derivation. 

Theorem 4.5. Given the general solution {(f)^ of the honeycomb linear system (|2.1ip - (|2.12p . and given its 
particular solution {0~^,9~), then the solution (0+,^^) of the system 



(4.16) 



(4.17) 0+ _ ^- ^ ^ {oZi.-n't'Zn ^I/iC/,-//; , 

(4.18) r - ^tjjj = ^ - OZj^Zj^^u) , 

with r given by equation (|3.2p . satisfies a new honeycomb linear system with the coefficients 

(4.19) a = -^eZnOZj^^n^ b = -^-^B-_,QZu-iv Sjjj = -^^,^1,,. 

Proof. The proof can be done by direct verification. However, in the spirit of the paper, we present its sublat- 
tice derivation. Consider the wave function ^ defined in the vertices of the hexahedron, like in Figure [B] The 
basic geometric property of the B-reduced fundamental transformation (see Remark 11.21) implies the possibility of 
eliminating the "black points of the hexahedron" from the discrete Moutard equation 

(4.20) 0-V-i-3 = /ii-3(V'-3-^^-i), 
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and the two Moutard transformation equations 

(4.21) = 

(4.22) ^_3_V,_i_3 = ^-i^(^_i_3-V'_3), 

c'-3 

which gives 

(4.23) 0^ - 0_i_3Vi-l-3 = /"i_3(^-3^-l - e-lV'-3). 

Define 0+ = 0_2_3V'-2-3, 0^ = ^-i-3V'-i-3 ill order to respect the duahty transformation, and notice that 
9^ ^ 9i, 6^ = 62 are solutions of the honeycomb linear problem (|2.1ip - (|2.12p . In such notation, and taking into 
consideration equation (|4.1ip . we can rewrite equation (|4.23p in the form of the first equation (|4.16p of the Darboux 
transformation system. A similar procedure gives equations 



-2-3^-2-3 - 6l_i_3Vi_i_3 = —pr^ (6'_3V'-l-2-3 - ^'-1-2-3^'- 



3 



-2-3 



^-1-3^-1-3 - 6'-l-2V'-l-2 = T23 (^'-1^-1-2-3 - 6'-l-2-3V'-l), 

/-1-2-3 

which, using equations (|4.12|) . give the remaining part (|4.17p and ()4.18|) of the Darboux tranformation system. 

Because tp and satisfy the same Moutard system, therefore {tpi,tp2) and (6*5"^,^^^) satisfy the same honey- 
comb afiine linear system. By Corollarv l2.41 the functions (0"*^, 4>~) satisfy the honeycomb affine linear system, the 
coefficients of which can be found from equations (|2.3p with in place of /'■' , supplemented by the (—1 — 2 — 3) 
shift and by the change (|2.16p with respect to the above-mentioned gauge with {a^,a^) = (^l2-3' ^-1-3)- This 
implies that the transformed coefficients read 



f-2-3' 


f-1-3 


(7-317-1- 


2-3 




-2-3 


fl2 

J-1-2 


-3 


9-1-2^ 


^-1-3 


O-16-I- 


2-3 




'2-3 


f23 

J -1-2^ 


-3 




3^^'-l-3 


= -f'-l 


-3^- 



SIJI 

which, using equations (|4.11|) and ()4.12p . take the form as in equation ()4.19p . □ 



5. Concluding remarks 

We have derived the integrability properties of the self-adjoint linear problems on the triangular and honeycomb 
lattices from the well known theory of the 3D system of discrete Moutard equations. This result can be the starting 
point for several research programs. We would like to point out some of them. 

In soliton theory, a linear problem admitting the Darboux type transformations can be used to construct evolution 
equations (with the time being discrete or continuous variable) compatible with the linear system; the standard 
example is the KdV hierarchy. As we have already remarked, the Laplace transformation equations (|3.1ip - (|3.12p 
provide such a nonlinear system. More complicated evolutions of the triangular and honeycomb lattices are being 
studied, see [55] , 

Given integrable reduction of the 3D Moutard system (for example, the generalized isothermic lattice [IS]), 
by restricting it to the quasiregular rhombic tiling one can obtain the corresponding integrable system on the 
triangular or on the honeycomb lattice. In particular, to obtain the constant angles circle patterns [S] one starts 
from further reduction of the generalized isothermic lattice. We remark, that such a process has already been built 
in the " consistency around the cube" property [SHI H] . 

Finally, using the "cut and project" method [35] in conjunction with the sublattice approach [21] one can derive 
linear problems (see also Remark 1 2. 3 1) and integrable systems on more complicated graphs. 
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